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ABSTRACT 

We compare the rotation rate of neutron stars in low-mass X-ray binaries (LMXBs) with the orbital 
period of the binaries. We find that, while short orbital period LMXBs span a range of neutron star 
rotation rates, all the long period LMXBs have fast rotators. We also find that the rotation rates are 
highest for the systems with the highest mean mass accretion rates, as can be expected if the accretion 
rate correlates with the orbital period. We show that these properties can be understood by a balance 
between spin-up due to accretion and spin-down due to gravitational radiation. Our scenario indicates 
that the gravitational radiation emitted by these systems may be detectable by future ground-based 
gravitational wave detectors. 

Subject headings: accretion, accretion disks — gravitational waves — pulsars: general — stars: 
evolution — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

Accreting millisecond X-ray pulsars are the evo- 
lutionary link between the (million-year old) clas- 
sical radio puls ars and the (billion-year old) mil- 
liseco n d pulsars dSmarr fc Blandford) Il976t lAlpar et al.l 
119821 Radhak rishnan fc Srinivasa r Hl98l . The ac- 
creting millisecond puls ars are found in low-m a ss X- 
ray binaries (LMXBs: IChakrabartv fc Morgan! 119981 : 
iWijnands fc van der Klisl 119981 ). systems in which the 
pulsar is accreting from a low-mass stellar companion. 
The details of how the transition from a rotating neu- 
tron star in a low-mass X-ray binar y to a millisecond ra - 
dio pulsar takes place are unknown (|Bhattacharvalll995l ) ■ 
although a system has been found rece ntly which shows 
evide nce for having made this transition ([Archibald et all 
[2001 . 

We show here that a significant clue regarding the tran- 
sition may be contained in the relation between the ro- 
tation rate and orbital period for the known accreting 
millisecond pulsars. This trend suggests that a common 
physical mechanism(s) may be responsible. We outline 
a simple scenario that can explain the evolution of the 
rotation rate and orbital period for accreting millisecond 
pulsars and show how this naturally leads to the observed 
population. The proposed evolutionary scenario predicts 
gravitational wave emission at a level that may be de- 
tected by next-generation ground-based observatories. 

2. DATA 

We use the tabulated values for the most rapidly- 
rotating (v s > 100 Hz) accreting millisecond pulsar 
(AMSP) spin frequency v s , orbital period P rb, and mass 
accret ion rate M from iGallowavl (|2008l ): iWatts et al.1 
([2008f FI. The spin frequencies of some AMSPs arc 
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derived fro m brightness osc i llations seen during X- 
ray b ursts ([Strohmaver et al.1 119961 : IChakrabartv et all 
l200l . which are due to thermonuclear burning of ac- 
creted material on the surf ace of the neutron star 
([Strohmaver fc BTld stcn 200(f). We do not consider spin 
frequencies derived from kHz quasi-periodic oscillations 
since the correlation between the spin an d quasi-periodic 
oscill ation frequencies is uncertain ([Wiinands et al.l 
l200l . 

The trend in the measured values of the spin frequency 
and orbital period for the AMSPs is shown in Fig.[U Also 
shown for comparison are values for the known (rotation- 
powered) millis econd radio pulsars from the ATNF Pul- 
sar Catalogu^l ([Manchester et al.ll2005[ ): these rotation- 
powered pulsars are older and thought to be the descen- 
dants of the AMSPs. While the rotation-powered pulsars 
are distributed fairly uniformly in the diagram of v s ver- 
sus P or b, the AMSPs are not. In particular, at the high- 
est rotation rates, the orbital periods of the AMSPs span 
the range from about 0.5 hr to 20 hr, but there are no 
AMSPs with 100 Hz < u s < 440 Hz and P orh > 4.3 hr 
A sim ilar plot to Fig. [TJ was presented by iKaaret et al.1 
((20061) . who remarked upon a possible absence of slowly- 
rotating, long orbital period systems; they suggested that 
this may be due to an accretion rate that varies with or- 
bital period but made no further discussion. A full pop- 
ulation synthesis of neutron star spin and orbital periods 
is beyond the scope of this work. Nevertheless, a quali- 
tative explanation of the ^ S -P rb distribution of AMSPs 
can be obtained by invoking simple physical processes 
that occur in their evolution. 

3. AMSP EVOLUTION 

In close binary systems such as the LMXBs, matter 
near the surface of the companion star is transferred to 
the primary star (the latter is a neutron star in the case 
of AMSPs) when the companion radius R c is equal to 
the Roche lobe size i?L- This size is the distance from 
the center of the donor star to the inner Lagrangian 
point, where the gravitational forces from the compan- 
ion and primary are equal and opposite, and is given by 

4 http:/ /www. atnf.csiro.au/research/pulsar/psrcat 
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Fig. 1. — Observed spin frequency as a function of orbital pe- 
riod for the accreting millisecond pulsars (triangles) and rotation- 
powered radio pulsars (stars). Th e spins and orbi t s for the ac- 
cretin g millisecond pulsars are from [Galloway (2008); Wat ts et al.l 
l|2008l ). whi le those for the radio pu lsars are from the ATNF Pulsar 
Catalogue (Manchester ct al. 2005). 



(|Paczvnskj|197H ) 

i? L ~ 0.46 a[m c /(m x + m c )] 1/3 , 



(1) 



where a is the orbital separation between the neutron 
star of mass m x and the companion star of mass m c 
(m x and m c are in units of solar mass M©); hereafter, 
we assume m x = 1.4. When the companion is a main- 
sequence star, m c ~ R c /R(t,. where Rq is the solar radius 
()Verbunt k van den Heuvelll 1995ft . Using Kepler's Third 
Law, 



O.8P (m x 



a/3 



(P orb /2 hr) 



2/3 



(2) 



we obtain a relationship betw een the mass of t he com- 
panion and the orbital period ([King et al.lll996|) . 



0.23m c (P orb /2 hr), 



(3) 



where m c = m c /m^ and m^ s is the mass of a main- 
sequence star that just fills its Roche lobe at P rb! 
the factor of m c (< 1) is due to the fact that the 
companion star in a binary has a lar ger radius for its 
mass than an isolate d star (see, e.g.. iSills et al.l [20001 : 
lAndronov et al.ll2003D and encapsulates our uncertainty 
in the (less evolved, since we only consider m c < m x ) 
evolutionary state of the companion. 

We now examine the various physical processes that 
would lead to a change in the angular momentum of the 
orbit or neutron star spin and thus cause a source to move 
within the ivP orb -planc of Fig. [TJ Mass transfer from 
the companion to the neutron star can cause an increase 
in the size of the orbit, while orbital angular momen- 
tum loss due to magnetic braking or gravitational wave 
emission causes the orbit to decrease. The orbital period 
which separates exp ansion and decay is estimated to be 
> 0.5 days (see, e.g.. lPod siadlows ki et al.l l2002: M a k Lil 



I2009D . For orbital periods longer than one hour, mag- 
netic braking is dominant. We consider two prescrip- 
tions fOT^_the_iriagnetic brakin g torque. The first one, 
from lVerbunt k Zwaan ( 198l]) , results in a timescale for 
decay 

r d v c L y « (2.1 x 10 8 yr) m x / 3 m- 4 / 3 (P orb /2 hr)" 2 / 3 (4 ) 

and a mass transfer rate M_n (in units of 
10- U M Q yr" 1 ) (see also King et "ail [T^96h 



MY& = 80m x 2 / 3 m c 7 / 3 (P orb /2 hr) 5 / 3 . 



(5) 



However, observations of rapidly rotating low-mass stars 
in open clusters suggest t hat the rate of an g ular m o- 
mentum l oss describ e d in iVerbunt k Zwaanl (1198 lT) is 
too high (ISills et al.l 120001: lAndronov et al.l 120031: see 
also lYungelson k Lasotal |2008[ for more discussion). 
Ther efore we also use th e magnetic brakin g torque 
from iChabover et all (Il995l): ISills et al.1 pOOOl ) (see also 
lAndronov et alj|200a IMa k Lill2009ft . which results in 

: (2.0 x 10 8 yr) m 2/3 TO c (P orb /2 hr) 7/3 (6) 
:87m x 2 / 3 (P orb /2hrr 4 / 3 . (7) 
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Note that the mass accretion rate given by cq. in- 
creases with orbital period, whereas that given by eq. ([7]) 
decreases. We also note that, at very short orbital pe- 
riods P orb ^ 4 hr, gravitational radiation becomes im- 
portant and some AMSP companions are degenerate 
stars. These systems may then be evolving to longer 
orbital periods, but they will evolve slowly because the 
mass transfer rate and angular momentum loss rate 
are low. For example, SAX J1808. 4-3658 has a de- 
generate companion and increasing orbital period with 
timescale P orb /Porb ~ 6 x 10 7 yr (IPods iadlowsk i et al.l 
12001 [DTSalvo et al.ll200i IHartman et al.ll2008l I2009D . 

Accretion of angular momentum-carrying matter 
from the companion can spin up the neutron 
star. The timescale for spin-up is r su ps (1.5 x 

10 9 yr)m x 3/7 P g " 2/7 AC 1 6 1 /7 (;y s /100 Hz), where B s is the 
neutron star magnetic field (in units of 10 s G). Using 
the mass accretion rate from eqs. © and ([7]), we obtain 

r s v u z = (3.4 x 10 7 yr) m l J 7 B^ 7 m-\P olh /2 hr)" 10 / 7 
x(^ s /100 Hz) 

= (3.2 x 10 7 yr) m x / 7 p- 2/7 (P orb /2 hr) 8 / 7 

x(i/ s /100Hz). (8) 

Spin-down of the neutron star by gravitational 
quadrupole radiation occurs on a timescale 

r sd « (2.9 x 10 10 yr) m x ^I^/lOO Hz)" 4 , (9) 

where e = 10~ 8 £_s is the quadrupole ellipticity. 

We ignore spin-down by electromagnetic dipole ra- 
diation, which only becomes dominant at > 
10 £_8(^s/100 Hz). It is instructive to estimate the mag- 
netic field that would be required for the AMSPs to be 
in magnetic spin equilibrium. Magnetic spin equilibrium 
occurs when the rotation period of the neutron star is 
equal to the Kcplerian period of th e inner accretion disk 
at the magnctosphcre boundary ((Davidson k Ostrikerl 
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Fig. 2. — Magnetic field of accreting millisecond pulsars obtained 
from assuming the neutron star rotates at spin equilibrium (see 
text) for the observed mass accretion rate. For a given M, neutron 
stars with B > £J cq are mass-propellers, while stars with B < B cq 
are accrctors. 



(10) 



Il973|k this leads to an equilibrium spin period 

^ s oq = (270 Hz) m 5 / 7 B~ 6/1 Aff£ 
or magnetic field 
P cq = (3.2 x 10 s G) m 5 J 6 M^^/lOO Hz)" 7 / 6 . (11) 



Neutron stars with v s > (or B > B cq ) can act as 
mass-propellers, as their rapid rotation creates a cen- 
trifugal barrier to accretion, while stars with v a < i/| q can 
accrete mass, gain angular momentum, and be spun-up. 

We consider whether the observed properties of neu- 
tron stars are consistent with the idea that they are 
in magnetic spin e quilbrium. First, AMSPs are likely 
to ha ve B < B eq ([Hartman et al.| [200l [201ll iPatrunol 
[20101: iPapitto et all 120 111) . Using the maximum accre- 
tion rate observed during outbursts (|Gallowavl I2008D . 
the high implied B cq (see Fig. [2]) could channel the ac- 
cretion flow and produce persistent coherent pulsations 
that are not detected; the absence of persistent pulsa- 
tions in the high accretion rate systems, where the spin 
is generally measured from burst oscillations, could be 
the result of magnetic scre ening by the accreted mate- 
rial |C umming et al. I [20011 ). In addition, some systems 
do not show changes to a propeller state at the accretion 
rates estimated from spin equilibrium dBarret fc Olivd 
l200l iMaccarone fc Coppil [200l [Gladstone et al.ll2007D . 

The evolution of v a and P orb is primarily determined 
by the process with the shortest timescale. A comparison 
of the timescales is shown in Figs. [3] and [4] Over the 
ranges displayed, there are three regions/regimes: spin- 
up from the accretion torque [sec cq. ([8])] , spin-down from 
gravitational quadrupolc radiation [see eq. ©], and orbit 
decay from magnetic braking [see eq. ((H) or ©]. An 
AMSP moves to higher v s (spin-up regime) if the pulsar 
rotation rate is below both v\ (from r su < r sc i) and v% 




Po* (hr) 

Fig. 3. — Evolution regimes of accreting millisecond pul- 
sars: AMSPs in the spin-up/spin-down regime (I/II) move to 
higher/lower spin frequencies, as well as to shorter orbital peri- 
ods but on a longer timescale; AMSPs in the orbit decay regime 
(III/IV) move to shorter orbital periods, as well as to lower/higher 
spin frequencies but on a longer timescale. The lines (labeled 
U 2 Z > U 3 Z ") separating the regimes are given by eqs. I12H - 114I . The 
open triangles denote AMSPs, while the solid triangles are the 
two AMSPs, SAX J1808.4-3658 at (P orb , u a )=(2 hr, 401 Hz) and 
IGR J00291+5934 at (2.5 hr, 599 Hz), which show short-term spin- 
up and long-term spin-down (see text). 



(from r su < T dccay ), where v x is given by 



,vz 



= (330 Hz)e_g / ° J B, 



-2/5 R 2/35.2/5 



(P orb /2 hr) 



2/7 



^ cs = (340 Hz) eZl /5 B 2 s /35 (P 01 . h /2 hr)- 8 / 35 
and v% is given by 

= (750 Hz) Bl' 7 rnll\P olh /2 hr) 16 / 21 
^ 3 CS = (740 Hz) B 8 2/7 m c (P orb /2 hr) 25 / 21 . 



(12) 



(13) 



If the rotation rate is above v\ and v-i (from T s a < Tdecay), 
where v-i is given by 



, 2 VZ = (270 Hz) e:^ /2 my 3 (P orb /2 hr) 1 ^ 



,cs . 



(280 Hz)e_ 



8 

-V2„ 



-1/4 



(P orb /2 hr) 



-7/12 



(14) 



then the AMSP moves to lower v s (spin-down regime). 
Finally, a decrease in orbital period (orbit decay regime) 
occurs when the rotation rate is below vi and above U3. 
We note that deviations from low valu es of rh c (up to ~ 
1, e.g., due to degenerate companions; King et all [19961) 
have the strongest effect on ^3, which do not change our 
results qualitatively. 

Though dependent on the various parameters, the ab- 
solute timescale of the dominant process in each regime 
bears out the observed population, with a slight pref- 
erence for a r su that decreases with P orb , such as rj^ z 
(which is due to a magnetic braking torque that scales 
with Pj rb , where 7 > 1/3). In the spin- up region (long 
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Fig. 4. — Evolution regimes of accreting millisecond pul- 
sars: AMSPs in the spin-up/spin-down regime (I/II) move to 
higher/lower spin frequencies, as well as to shorter orbital peri- 
ods but on a longer timescale; AMSPs in the orbit decay regime 
(III/IV) move to shorter orbital periods, as well as to lower/higher 
spin frequencies but on a longer timescale. The lines (labeled i/f s , 
j/^ s , v^) separating the regimes are given by eqs. H12H - II14I I, The 
open triangles denote AMSPs, while the solid triangles are the 
two AMSPs, SAX J1808.4-3658 at (P orb , u s )=(2 hr, 401 Hz) and 
IGR J00291+5934 at (2.5 hr, 599 Hz), which show short-term spin- 
up and long-term spin-down (see text). 

Porb, low V s ), 

t™~(7 x 10 6 yr)m c - 2 (F orb /10 hr)- 10 / 7 (^ s /200 Hz) 

rg s ~(4x 10 8 yr)(P orb /10 hr) 8/ > s /200 Hz). (15) 

In the spin-down region (high z/ s ), evolution is indepen- 
dent of orbital period and occurs on a timescale 



r sd ~(2x 10 9 yr)( £ /10- 9 )- 2 (^ s /500 Hz)" 



(16) 



In the orbit decay region (short P or b), evolution is inde- 
pendent of spin frequency and occurs on a timescale 

T d v J ay ~(2 x 10 8 yr)m c - 4 / 3 (P orb /3 hr)- 2 / 3 

7d C ecay ~(6x 10 8 yr)m c (P orb /3 hr) 7 / 3 . (17) 

The two AMSPs, SAX J1808.4- 3658 (IHartman et al l 
20091) and IGR J00291+5934 jlFalanga et all 120051 
Patruno|[2010t IHartman et~aT]|20 lit IPapitto et al.ll2011[ ). 



show an overall spin-down with timescales f s / 1 | ~ 
2 x 10 10 yr and 5 x 10 9 yr, respectively (interupted by 
short spin-ups during outbursts with timescales 10 s yr 
and 2 x 10 7 yr, respectively). 

4. DISCUSSION 

We can now understand the observed population of 
AMSPs (especially the absence of sources at low spin 
and long orbital period) as a result of the evolution of 
v s and P or b, depending on the process (magnetic brak- 
ing, mass accretion, and gravitational radiation) with the 
shortest timescale. AMSPs born at low v s and long P or b 
very quickly spin up to high v s . Then on much longer 



timescales, these fast spinning sources slow down, and 
their orbits decrease. Once mass accretion ceases, there 
is no longer a spin-up torque; the binary then contains a 
rotation-powered millisecond pulsar that can move into 
the (former) spin-up region in Figs. [3] or [4] by spinning 
down (as a result of gravitational wave or electromag- 
netic radiation) or expanding its orbit (see Fig. [T|). 

We utilize two very different prescriptions for angu- 
lar momentum loss due to magnetic braking and find 
that we can vacate the low spin, long orbital period 
spin-up region for reasonable parameter choices in both 
cases. We note a slight preference for the results ob- 
tained using eqs. (@| and (JS|), as compared to eqs. (O 
and flTJ). Future observations, for example, measurement 
of the mass accretion rate dependence on orbital period, 
could distinguish between the two models. There are 
presently six AMSPs with unknown P or b, three of which 
have v s < 363 Hz. Measurement of a long orbital pe- 
riod for these systems or discovery of new systems at 
the relatively low v s (100 Hz < v s < 400 Hz) and long 
Porb would challenge the scenario proposed here. Clearly 
more detailed studies of short orbital period systems 
are required to address the trend seen, as well as ac- 
counting f or spin evolution in binary evolution models 
(see, e.g., IPodsiadlowski et all 12002 1: ILamb fc Yul 120051: 
iFerrario fe Wickramasinghell2007f r 

The most rapidl y-rotating millisecond radio pulsar 
(|Hessels et all [20061 ) and AMSP have v s = 716 Hz and 
620 Hz, respectively. These frequencies are far below 
the theoretical maximum (at > 1 kHz), a bove which 
the centrifugal force caus es mass-shedding (|Cook et all 
119941: lHaensel et all 1 1 9991 ). It is thought that angular 
momentum loss from gravitational radiation could be re- 
sponsib le for a spin limit below the mass-shedding max- 
imum (|Wagoner| [T984: Bil dstenl 119981 : lAndersson et all 
119991: IMelatos fe PavneH2005D . However, theoretical pre- 
dictions for this limit are very unce rtain, with only the 
current gravitational wa ve searches (lAbbott et alll2010D 
and X-ray observations (|Chakrabartv et all 120031 ) serv- 
ing as constraints. In our scenario, the observed posi- 
tion of AMSPs in Figs. [3] and H relative to the differ- 
ent evolution regions is suggestive that the amplitude 
of the mass quadrupole that produces gravitational ra- 
diation is e > 10~ 9 . At lower e, spin-down by gravi- 
tational radiation becomes irrelevant [see eqs. (fT2j) and 
(|14p], and all AMSPs would be in spin-up or orbit expan- 
sion or decay, contrary to what is seen (IHartman et all 
[2001 [20TT1 : IPatrunol[2010l: IPapitto et alll2011D : note that 
spin-down by electromagnetic radiation gives a similar 
region in Figs. [3] and @] only if B > 10 9 G. Our estimated 
quadr upole ellipticity is far belo w the theoretical maxi- 
mum ([Horowitz fe Kadaul[2009T) of 4 x 10" 4 and below 
the curren t limit of 7 x 10~ 8 set by gravitational wave 
detectors ([Abbott et al.ll2010D. Though observationally 
challenging (see. e.g. JWatts et alll2008[) . future searches 
by Advanced LIGO o r the proposed Einstei n Telescope 
([Punturo et al.N2010t lAndersson et~aTl 120111) could pro- 
vide direct evidence for the evolutionary scenario out- 
lined here. 
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